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The coherent states are reviewed with particular application to the free particle 
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CT 1 . I. INTRODUCTION 

Coherent states were discovered by Schrodingerjdi] in 1926 but remained unex- 
^ \ ploited until Glauber [2j rediscovered them in 1963 and developed a profound for- 

■ malism for them with important applications; for instance, the "Glauber states" are 

fundamental in the description of the electromagnetic field for quantum optics. He 

O 

was awarded the Nobel price in 2005 for the development of the quantum theory 
of optical coherence. Besides its important applications in quantum optics, theses 
states provide an elegant mathematical description of the harmonic oscillator and 
in many textbooks they are presented as a didactic tool. However, the didactic ad- 
vantage of the coherent states has been ignored in the treatment of a free particle. 
It is, in fact, very convenient to present the gaussian states of a free particle, not 
as is usually done by assuming a gaussian state in momentum space and (through 
Fourier transformation) obtaining the corresponding gaussian in configuration space 
with a quite unappealing mathematics, but instead, one can present the coherent 
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states for a free particle, and with simple and elegant algebraic methods, derive all 
properties of the gaussian states. The treatment of the coherent states of a free par- 
ticle is similar to the corresponding treatment of the harmonic oscillator but there 
are remarkable differences and therefore it is convenient to present the complete 
derivation that can be directly used in the teaching of quantum mechanics in the 
advanced undergraduate or graduate level. 



II. THE FREE PARTICLE SYSTEM 



Let us consider the simplest quantum system consisting in one free structure- 
less particle moving in a one dimensional space. Let H be the Hilbert space for 
the description of the system. It is usual to choose the space of square integrable 
functions £2 for this space but in principle it can remain abstract. Position X and 
momentum P are the unique independent observables of this system. Their eigen- 
vectors {fx} and {4> p }, build two mutually unbiased bases and their commutation 
relation is [X, P] = ih. (We assume here an infinite dimension of the Hilbert space; 
for finite dimension, in a lattice for instance, the commutation relation is more 
complicated^]). The hamiltonian of the system is H = J^- and the time evolution 
operator is given by U t = exp(—iHt). 



III. THE COHERENT STATES 



It is very convenient to introduce a length scale A for the system that will allow 
us, together with h, to make position and momentum dimensionless in order to be 
able to add them as simple numbers. For this scale we can choose whatever we like, 
for instance the Compton length of the particle or any measure of derealization of 
the particle in space (we will later see that it is convenient to choose A proportional 
to the width of the space distribution). With this scale we define the operator A 
and its hermitian adjoint 

A ~ Avl X ~ l hV^ P ' ^ 
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These operators are dimensionless and their commutation relation is 

[A, ^] = I . (3) 
The relations above can be trivially inverted to obtain 

Now we can define the coherent states ip a as the eigenvectors of A: 

Aip a = ai\) a . (6) 

Since A is not hermitian its eigenvalues are complex: a G C. It is remarkable that 
we can very easily calculate these eigenvalues, related with observable properties of 
an ensemble of particles in this coherent state. Let us find (X), (P), (X 2 ), (P 2 ) and 
the widths A x and A p for the particle in a coherent state ip a - 

(7) 

and similarly we obtain 

<P> = -I; («-*•). (8) 

Adding them we obtain the eigenvalue a given in terms of the expectation values of 
position and momentum 

We could have obtained this relation simply by taking the expectation values in 
Eq.(l). The coherent states are then completely determined by the expectation 
values of position and momentum. For (X 2 ) we have 

(X 2 ) = —((A + At) 2 ) = — (A 2 + At 2 + Atf + A^A) 

= ^ {A 2 + At 2 + 2AU + 1) = ^ (a 2 + a* 2 + 2a*a + l) , (10) 
and similarly we get 

(P 2 ) = ~ {a 2 + a* 2 - 2a*a + l) . (11) 
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The widths of the position and momentum distributions are then 

A* = (X 2 ) - (X) 2 = y , (12) 

K = (p 2 ) - (P) 2 = ^ , (13) 

where we see that the length scale is related to the width in position and we get 
the important result that, in a coherent state, the uncertainty relation is optimized 
A x A p = §. For this reason, it is often mentioned that the coherent states are the 
quantum states closest to classical behaviour. In these coherent states, position 
and momentum have the least correlation (the real part of the quantum covariance 
function^ vanishes) in the sense that (XP + PX) - 2(X)(P) = 0. 

We have seen that the eigenvectors ip a of the operator A have very interesting 
properties and we may wonder whether the eigenvectors of A* may be of some 
interest. We can see however that the eigenvalue equation A^xp — 0Xp does not 
have solutions because, if xp would exist then we would obtain an absurd result. In 
fact, if we assume that xp exists, applying the same techniques as before, we can 
calculate (X), (X 2 ) and the width A 2 = (X 2 ) — (X) 2 for the particle in the state xp- 
Doing this we obtain A 2 = — A 2 /2; a negative result! This is of course impossible 
because A 2 is the expectation value of the positive operator (X — (X)) 2 , therefore 
there is no xp satisfying A^xp = Pxp- 



IV. TIME EVOLUTION OF COHERENT STATES 

So far, the treatment of the coherent states for a free particle and for the har- 
monic oscillator are similar and we can now present some important differences. In 
the harmonic oscillator, one defines a "number operator" X = AA* having non- 
negative integer eigenvalues, with A and A* as shift operators for the eigenvectors 
of X. This operator is very important for the harmonic oscillator because it has a 
relevant physical meaning; it is, essentially, the hamiltonian operator that controls 
the time evolution. Of course we can also define this number operator for the free 
particle, and we will do so in the next section, but in this case it is not related 
to a physically interesting observable. In the free particle case, the hamiltonian is 
given by the momentum squared and the energy eigenvectors are the same as the 
momentum eigenvectors (f) p (notice that the energy eigenvalue E = p 2 /2m has a 
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twofold degeneracy with eigenvectors P and <p- p .) 

If at time t = the free particle is in the coherent state ip a then at time t it will 
be in the state given by the time evolution operator 

*(t) = exp {-'-Ht) *. = exp glf) «. = exp - . 

(14) 

Here we can see another important difference with the harmonic oscillator. One can 
prove that the coherent states of the harmonic oscillator remain coherent through 
the time evolution whereas in the free particle case the time evolution destroys the 
coherent states. That is, ^(t), for t ^ 0, is no longer a coherent state. This is phys- 
ically understood because, as is well known, the width in the position distribution, 
A x , increases with time for the free particle and the width in momentum, A p , re- 
mains constant (because the momentum distribution is time independent); therefore 
the uncertainty product A^Ap will also increase with time and will no longer have 
the minimal value h/2 characteristic of the coherent states. Anyway we will now 
prove that \I/(t) is not an eigenvector of A and therefore it is not a coherent state. 
For this, consider the time evolution applied to the Eq.(j6]) that results in 



exp ( -~Ht \ A exp f -Htj V(t) = atf(t) . (15) 

Now, for the operator in the left hand side of this equation we can use the general 
expression valid for any constant 7 and any operators A and B 

exp (- 7 5) A exp ( 7 fl) = A+^B, A]+^[B, [B, A}}+- ■ -+^[B, [B, A}} •••]+•• • 

(16) 

To prove this equation, consider the function /( 7 ) = exp (— jB) A exp (7-B) and 
calculate the derivatives with respect to 7. Now making a Taylor series expansion 
about 7 = we obtain the relation above. In our case, however, the series is 
interrupted after the second term: 

[H,A] = — [P 2 ,A] = — I — [P 2 ,X] = ^-P, (17) 



and [H, [H,A]] = 0. Therefore Eq. (jl5p becomes 

t 

V2m\ 



A + -jj — P ) *(t) = aW(t) 
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that is, 

A* (t) = — Z—PV(t) + aV(t) ^ pV(t) . (19) 
V2mA 

The last inequality follows because ^(t) is not an eigenvector of P (neither is ip a ) 
and therefore it is also not an eigenvector of A. 



V. THE NUMBER OPERATOR 



As was announced at the beginning of last section, we will now study the number 
operator N = AA^ that has interesting formal properties although it is not related to 
any important observable of the free particle physical system, whose most relevant 
observables are position and momentum. Another operator similar to N is A' A and 
one may wonder whether it may bring something interesting. However, due to the 
commutation relation [A, A'] = I, it is clear that this operator is not significantly 
different from N and it amounts only to the addition of a constant. The operator 
N is hermitian and positive N > 1, that is, it is bounded from bellow. This is 
important because one can prove that for such an operator there exist eigenvalues 
and eigenvectors (A n ,x„). The prove that = N is trivial and the positivity 
follows from 

(\I>, AA^) = (A^,A^) = \\A^\\ 2 > ,Vtf G U . (20) 

Then we have Nx n = ^nXn, with \ n > 0. Using the commutation relations 

[N, A^] = A ] and [N, A] = -A (21) 

one can easily prove that A^Xn is another eigenvector of N but corresponding to 
the eigenvalue A n + 1 and A\ n is another eigenvector of N but corresponding to the 
eigenvalue A n — 1. Now, if we start with an arbitrary eigenvector Xn and apply A a 
sufficient large number of times, then we would get an eigenvector with a negative 
eigenvalue in contradiction with the positivity of N. The solution of this difficulty is 
that the eigenvalues must be integer numbers. Then we can set A n = n = 0, 1, 2, . . ., 
and there is a lowest stair Xo such that Axo = from which we can obtain any Xn 
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by applying n times the raising operator A'. Summarizing, we have: 

Nxn = n Xn with n = 0, 1, 2, . . . (22) 

A^xn = Vn + 1 Xn+i , (23) 

Axn = Vn Xn-i , in particular, Axo = , (24) 

Xn = -^(A^Txo • (25) 



Notice that from Eqs.(j6]and 24) it follows that the lowest number state coincides 
with the coherent state for a = 0, that is, Xo = V'o- We will later see how every 
coherent state i/j a can be expanded in the number states basis {Xn}- 

VI. THE DRIFT OPERATOR 

An interesting operator for the coherent states formalism is the drift operator 
defined as 

D(a) = exp (aA ] - a* A) 

= exp ^— ^ij - ^ ex P ( a ^) exp (—a* A) 



exp ^"^j - ^ ex P (~ ®* A) exp (aA^) . 



(26) 



The last two expressions for D(a) follow from the Hausdorff-Baker-Campbell (HBC) 
relation: 

exp(i?) exp(S) = exp (R + S) exp (^y^J > valid if ft i R > S ]] = t 5 ' = • 

(27) 

Repeated application of this identity results in 

exp(i?) exp(S) = exp(S) exp(R) exp([R, S]) , valid if [R, [R, S}} = [S, [R, S}} = . 

(28) 

With the HBC relation, one can also prove that the drift operator is unitary, 
Dt(a)D(a) = D(a)D^(a) = I, therefore D\a) = D _1 (a) = D(-a). The set 
{ Va 6 C} is a group whose product can also be calculated with the HBC 
relation as 

D{a)D{(3) = exp (^Lz^£\ D(a + (3) . (29) 
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We will prove that the operator D(a) can be used to relate coherent states corre- 
sponding to different values of a. More precisely, if ip is the coherent state for a = 
then we have 

i> a = D(a)?p . (30) 

In order to prove this, we will need some algebraic properties of the operators A and 
A\ By mathematical induction one can prove that, from the commutation relation 
[A, At] = I, it follows that [A, A^ n ] = nA^ n ~ l \ and from this, we can show that for 
every function F that can be expanded as a power series it is 

[A,F(Ai)] = ?*yp-. (31) 

Applying this to the drift operator we have 

[A, D(a)] = aD(a) . (32) 

With this result we can easily see that D(a)ipQ is an eigenvector of A corresponding 
to the eigenvalue a. In fact, we have: 

AD{a)ip = D{a)Atp + aD{a)tp = aD(a)^J . (33) 

Therefore D(a)ipo — kip Q , where k is some proportionality constant. In order to 
determine it, consider the norm of the equation 

\k\ 2 (^ a ) = (D( a )^ ,D(a)^p ) = {ip ,D\a)D(a)^ ) = (^o) = ||^o|| 2 • (34) 

Assuming that the coherent states are normalized, (i> a ,^ a ) — ll^all 2 — 1> then 
\k\ 2 = 1 and D(a)ipo and ip a can differ by a constant phase that can be set equal 
to 1 and therefore we obtain Eq. (l30l . Notice that the algebraic proof given here is 
different from the usual proof that involves the expansion of the coherent states in 
the basis of the eigenvectors of the number operator. In this approach we don't use 
such a basis because they are physically uninteresting for the free particle system 
and we prefer to rely only in the algebraic structure of the operators. The drift 
operator has an intuitive meaning that becomes evident when we write D(a) in 
terms of the position and momentum operators: Using Eqs.(l, 2, 9, 26), we have 

D(a) = e*p(j:({P)X-{X)P) 

= ^ (-i^) exp [n {p)x ) exp [-i {x)p 

= exp (±_(X)(P}\ exp (-I(X)P) exp (j:<P)x\ , (35) 



9 



and recalling that exp (— iaP) is a translation operator that performs X — > X + a 
and exp ^gX) changes the momentum P — )■ P + g, we see that D(a) transforms 
a coherent state corresponding to a particle at rest, (P) = 0, located at the origin, 
(X) = 0, into a coherent state for the particle at position (X) and moving with 
momentum (P). Notice that the different choice of what transformation is performed 
first, amounts only to a different constant phase. It is then interesting to notice that, 
effectively, D(a) performs a Galilei transformation on the state i/jq and therefore the 
group build with the drift operators D(a) is isomorph with the Galilei group. In 
fact, the phase appearing in Eq. fj29|) shows that this group is a projective group (not 
simple) as is the case with the Galilei group. The emergence of the Galilei group is 
an additional "bonus" of the formalism of coherent states. 

The arguments presented in this section are valid for a free particle as well as for 
the harmonic oscillator, or for any other potential. The only difference is that, in 
the harmonic oscillator case, the state ipo has the unique feature of being the ground 
state of the Hamiltonian and one may prefer to choose this state in order to generate 
all other coherent states. In the free particle case, ipo is a state as good as any other 
because of the Galilei invariance of the free particle system and one might prefer to 
write Eq. fl30l) in an unbiased way as ipp = D((3 — a)ip a . 



VII. COHERENT STATES ARE GAUSSIAN 



One of the most beautiful features of the Hilbert space formalism of quantum 
mechanics is that almost all results can be obtained without specifying a particular 
representation of the Hilbert space. Everything follows from the algebraic and geo- 
metric structure of the Hilbert space. Some features become however easier to grasp 
in some particular representation of the Hilbert space. For instance if we are inter- 
ested in the position distribution of a particle it becomes natural to choose the space 
£2 of square integrable functions of a real variable x. As is well known, the position 
operator X in this space is the multiplication by x and the momentum operator P 
is the derivative —^4z- In order to find the coherent states ip a (x) = ((p x , ip a ) in this 
space we must write and solve Eq.flH]) in £ 2 - That is 

ip a {x) = aip a (x) . (36) 



— -= x + i — -= i—ih—) 
\V2 hV2 dx' 
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In terms of a dimensionless variable y = x/(a/2A) this becomes 




(37) 



that is, 



dip, 



a 



2(a-y)dy . 



(38) 



Ipa 



Integrating we have 



ln(^ a ) 



y 2 + lay + Const. , 



(39) 



that is, 



ip a (y) = Cexp (-y 2 + 2ay) 



(40) 



The next step is to complete the square in the exponential and absorb the corre- 
sponding factor in the constant C (that will be anyway determined by normaliza- 
tion). After replacing the original variable x and a given by Eq.® with (X) = xq 



We see therefore that the coherent states of a free particle are the gaussian states (in 
order to appreciate the didactic advantages of the coherent states for a free particle, 
compare this derivation with the presentation of the of gaussian wave packets as is 
usually done in textbooks). 

One last remark in this section is that the coherent states in momentum repre- 
sentation, that is, the functions ip a (p) = ((p P ,ipa) are also gaussian. There is in fact 
a symmetry in the exchange of position and momentum for the coherent states. 

VIII. EXPANSION IN NUMBER STATES AND IN COHERENT STATES 

A very important feature of the Hilbert space formalism is the fact that the 
eigenvectors of an hermitian operator, that may be related to a physical observable, 
build a basis. If we expand the state ^ of a system in this basis, we may interpret 
the modulus squared of the expansion coefficients as the probabilities associated 
with the eigenvalues of the observable. 

A relevant question is then, if we can use the set of coherent states {ipa} in 
order to make an expansion of a state \l/ and interpret the expansion coefficients 



and (P) = po, and replacing the length scale A = \/2A x we obtain 



M*)= [27rA;^ 1/4 exp( 




(41) 
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\(ip a , ^)| 2 as some probability distribution. The first remark in order to answer 
this, is that the coherent states are the eigenvectors of an operator that is not 
hermitian and they do not build an orthonormal basis. They are not even linearly 
independent (if they were independent we could transform them into a basis by 
Schmidt orthogonalization procedure). We will show that the coherent states are 
not orthogonal, that is, (i[) a ,il)p) ^ 0. 

(Vo,^) = P(a)Vo,-D(^)^o) • (42) 
Using Eq.(26) and considering that exp(aA)-?/>o = ipo we get 

i^aitpp) = /exp f - ^rl ex P { aA] ) V>o, ex P exp (^) 

= exp f-^f - I^L) <^ , exp (a* A) exp (fitf) ^ > • (43) 



Now, using Eq. fl28|) we can permute the operators and we get 

2 ^ + j (i/>o, exp (f3A j ) exp (a* A) ip ) 

= exp I — + a P 

The three operators X, P, N define tree bases of the Hilbert space {<p x }, 
and {Xn}- The first two are continuous x,p G R and unbiased, that is, 1(^,^)1 is 
a constant independent of a; and p. The third basis is numerable and biased with 
the other two. One can calculate, for instance, that (<p x , Xn) = Xn{%) is related with 
the Hermite polynomials. We will not reproduce here this calculation that can be 
found in many quantum mechanics textbooks. Instead of this, we are interested in 
the expansion of the coherent states in these bases. We have already seen in section 
VII that the Fourier coefficients of the expansion of the coherent states in the bases 
of position and momentum, that is ((p x ,ip a ) and (<f) p , ip a ) , are gaussian functions. 
The corresponding coefficients, (Xm^a), of the expansion of the coherent states in 
this new basis {Xn} are given by 

(XnM = (-±=(Ai) n xoM = -^=(xo,A n ip a ) = 
vn! vn! 



1 a n a n ( |ct| 2 \ 

(^o,« n ^) = ^=(^o,^) = ^=exp -■ • ( 44 ) 



nl ynl \/n 
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We have therefore 

2 



^ a = Yl ~ri exp — =r ) Xn ■ ( 45 ) 



We can now go back to our question of whether it is possible to make an expansion 
of any state in terms of the coherent states. We have seen that these states are not 
orthogonal; however they build an (over) complete set of states in the sense that 

7T- 1 [d 2 a Va(VW> =1 (46) 

(we use here the correct Hilbert space notation for a projector: •); for those 

addict to the Dirac notation this is | 1 J r )( 1 i r |) and we can use this relation in order to 
expand any element \I/ as 

$ = I* = tt-1 f d 2 a i/} a (i/i a , *) . (47) 
./,,-. . 

The proof of the completeness relation in Eq. (1461) follows from the completeness of 
the basis {x-n}'- 

f f a n ( \ot 2 \\ a* m ( \ot 2 \\ 

\ d 2 a ip a (ip a , ■) = / d 2 a V -= exp — V -= exp — Xn(Xm, •) 

JadC JadC „ Vn! V Z J m V ™ \ 2 / 

= 5"xn(x m r)^=^ / d 2 a a n a* m exp (-|« 2 |) 
Vn!m! 

= 7T^Xn(Xn,-)=VTl. (48) 

n 

The integral in the complex plane is performed in polar representation a = r exp(id). 

The completeness relation in Eq. fl46|) can be used to expand any Hilbert space 
element. For instance, we can invert Eq. (1451) and give the number states expanded 
in the coherent states, using Eq. (l44"|) : 



Xn = EXn = 1 / d 2 a ^= exp ( ) ip a . (49) 



This expression, when written in the Hilbert space of the squared integrable func- 
tions, £2, allows an interesting representation of the Hermite polynomials in terms 
of gaussian functions. 

One last comment concerning the expansion in terms of the coherent states is 
that the set of elements {ip a }, as well as the basis elements {Xn}, are well defined 
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in the Hilbert space. Strictly speaking, the bases associated with position and 
momentum {tp x } and {<fi p } do not belong to the Hilbert space. In fact, if we choose 
the Hilbert space of the square integrable functions of some real variable, y, in the 
position representation where <p x (y) = S(x — y) and <p p (y) = -^exp^py), we can 
immediately see that they do not belong to the Hilbert space because they are not 
squared integrable. For this reason, in advanced quantum mechanics books, the 
Hilbert space is extended to the rigged Hilbert space that includes these elements 
that may be used as bases. This is precisely what we do when we calculate a Fourier 



transformation. The mathematics to justify the extension of the Hi 
the rigged Hilbert space is well explained in Chapter 1 of reference 



bert space to 



6]. 



IX. THE HUSIMI DISTRIBUTION 



The main interest in the expansion of a state \I/ of a system in terms of the eigen- 
vector of some observables is that the modulus squared of the expansion coefficient 
have the physical interpretation as the probability distribution of the correspond- 
ing eigenvalues. So, p{x) = \(ip x , \I/)| 2 and w(p) = |(0 P , \I/)| 2 are understood as the 
probability distributions for position and momentum (although this is perhaps a 
missnommer fjj ) . Similarly, we may try to interpret |(^q,,\I/)| 2 as some probability 
distribution for a. From the mathematical point of view, this seems to be possi- 
ble because, using the over-completeness relation Eq. (|46|) one can prove that the 
expectation value of any function F(A) can be given as 

(V,F(A)V) = (*,F(A)W) = ir~ 1 f d 2 aF(a) |(^ Q ,^)| 2 , (50) 

and therefore ir~ 1 \(%l) on \1/)| 2 can be interpreted as the probability for the realization 
of the value a. However, from the physical point of view this does not seem to make 
sense. The reason for this is that A is not an hermitian operator and therefore the 
complex eigenvalues a can not be associated with the result of a single measurement 
of some observable. Consequently, it is meaningless to say that the system has some 
value of a with some probability. In other words, let us consider the complex 
variable a as a function of two real variables x and p as in Eq.([9]) a(x,p) = -^^x + 
i^nP- Clearly, x and p can not be the result of any measurement of position and 
momentum because these can not be measured simultaneously. Neither is there an 
observable, with its corresponding hermitian operator, whose single measurement 
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provides simultaneously the expectation values of position and momentum. These 
expectation values can only be obtained in a multiplicity of measurements in an 
ensemble of identically prepared systems but not in some single measurement of one 
system. More precisely, we can see that given an arbitrary state that is, any 
Hilbert space element, with position expectation value (X) = then there 

exists no hermitian operator X such that X^> = (X)ty. Clearly this is impossible 
because if such an operator would exist, then all Hilbert space elements would belong 
to the basis associated with X. 

Even though \(i/i a , ^)| 2 is not a probability density, this quantity has been stud- 
ied in one of the several attempts to establish a compound probability distribution 
for position and momentum (see for instance Chapter 15 of reference jsj). That is, 
some function of x and p such that when integrated over p results in the position 
distribution p(x) and when integrated over x delivers the momentum distribution 
vj(p). The most famous attempt is the Wigner distribution that has these two prop- 
erties but is not everywhere nonnegative. The function pn(x,p) = 7r _1 | (i/j a (x,p), ^)| 2 
is known as the Husimi distribution that is everywhere nonnegative but, when in- 
tegrated over one of its variables, does not provide the correct distribution for the 
other one. There are interesting properties for this phase space "distribution" that 
can be found in ref . jfj and a profound comparative study and applications of many 
phase space distribution functions is found in ref.0]. 



X. CONCLUSIONS 



We have seen that the coherent states for a free particle, very well known for the 
harmonic oscillator case and in quantum optics, provide a very useful didactic tool 
that has not been exploited in quantum mechanics textbooks. For instance, the very 
elegant formal treatment discovered by Glauber allows a much nicer presentation of 
the gaussian wave packets for a free particle. The similarities and differences ap- 
pearing in the harmonic oscillator and in the free particles are emphasized providing 
a global outlook of quantum mechanics. Another advantage of the formalism is that 
one can easily relate the drift operators with the Galilei group and see that it must 
be a projective group. The expansion in terms of Hilbert space bases related with 
observables of position, momentum and number are exposed and compared with the 
expansion in the over-complete set of coherent states that provide a natural way for 
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introducing the Husimi distribution. The paper is written with complete calcula- 
tions, or with sufficient hints for the calculations, in a way that it can be directly 
handed to students as complement of any textbook. 
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